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ON THE RELATIONSHIP BETWEEN ACOUSTIC ENERGY DENSITY FLUX 
NEAR THE JET AXIS AND FAR- FIELD ACOUSTIC INTENSITY 


By Lucio Maestrello 
Langley Research Center 

SUMMARY 

By measurement and analysis, the relationship between the distribution of the out- 
flow of acoustic energy over the jet boundary and the far-field intensity is considered. 

The physical quantity used is the gradient of the pressure evaluated on a geometrical 
plane at the smallest possible radial distance from the jet axis, but outside the vortical 
region, in the area where the homogeneous wave equation is reasonably well satisfied. 

The numerical and experimental procedures involved have been checked out by using a 
known source. Results indicate that the acoustic power output per unit length of the jet, 
in the region from which the sound emanates, peaks at approximately 9 diameters down- 
stream. The acoustic emission for a jet Strouhal number of about 0.3 exceeds the emis- 
sion for all other Strouhal numbers nearly everywhere along the measurement plane. 
However, the far-field peak intensity distribution obtained from the contribution of each 
station was found to depend on the spatial extent of the region where sound emanates from 
the jet, which, in turn, depends more on the far-field angle than on the Strouhal number. 
For small far-field angles, from the jet axis, the frequency spectrum peak is independent 
of the Strouhal number, whereas at angles near 90° from the jet axis, it is a weak function 
of the Strouhal number. These results can also provide guidance on sound suppression by 
means of surface shielding which can only be achieved when the shielding plane extends 
beyond 20 nozzle diameters. 


INTRODUCTION 

The problem of jet noise must be formulated in a fashion amenable to experimenta- 
tion. The direct approach is the evaluation of Lighthill's integral (ref. 1) where the 
acoustic and the turbulent parts do not appear to be readily separable. During the past 
20 years, attempts have been made to evaluate this integral, on the basis of obtaining 
acoustic far-field information, but it has proved to be a formidable task. An alternative 
approach, suggested by Maestrello and McDaid (refs. 2 and 3), was to infer as much as 
possible about the sound generation region of. the jet from measurements made outside. 
This approach has the advantage of simplifying the experimental procedures; however, the 



simplification is achieved only at the expense of some loss of information. Lilley et al. 
(ref. 4) follow a theoretical approach similar to that of references 2 and 3 for the sound 
field radiated from a jet; in addition, they develop a theory for a cyclindrical plug-type 
flow for the sources in the jet, which has been interpreted as a special case of Doak's 
unified theory. (See ref. 5.) 

The chief objective of the present investigation is to determine the distribution of 
^ the acoustic energy density flux out of the jet in terms of the cross correlation of the 
pressure gradient measured on a geometrical plane outside the vortical region of the jet 
by following the approach of references 2 and 3. The acoustic energy density flux, 
although it does not specify exactly the sound generation region of the jet, does determine 
the zones from which the sound emanates. Furthermore, from the correlation of the 
pressure gradient on the plane, the far -field acoustic intensity is obtained, in a form 
which can be seen to be made up of contributions from various zones; thus, the gap 
between the near -field sound and the far -field sound is bridged. 

The approach of this paper presented two potential difficulties: (1) how close can 
one approach the jet and still remain in the acoustic region? and (2) what is the influence 
of the entrained flow on the energy flux ? To answer these questions, it was necessary to 
map the velocity field surrounding the jet and to establish the region in which the velocity 
fluctuations satisfy the homogeneous wave equation. 

This approach also permits checking the total acoustic power and the acoustic 
directivity since both are evaluated by the superposition of the sound field from the vari- 
ous stations along the jet. This verification of the accuracy of the method and of the 
results is extremely important especially since the method involves making, in effect, 
radiation field measurements relatively close to the jet; the building of a detailed picture 
of the sound field as it exits from the jet would thus be permitted, and also the manner in 
which it is redistributed into the far field can be determined. 

In the work reported here, the sound field emanating from an axisymmetric jet has 
been related to the far -field intensity. In addition, from the measurement of the cross - 
correlation function of the pressure gradient on the near plane, it has been possible to 
identify some of the modal features of the sound field emanating from a certain, well- 
defined conical surface enclosing the jet, this conical surface being defined by wrapping 
the measurement plane around the jet. This work may help to clarify some of the results 
obtained in the well-known experimental investigations by Mollo-Christensen (ref. 6) and 
by Crow and Champagne (ref. 7) who described a wavelike structure of the jet turbulence 
containing components of all frequencies, the different frequency components preserving 
their phase relationships over a few diameters. It may also quantify the more recent 
analytical work of Hardin (ref. 8) in which the source distribution produced by the orderly 
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structure of the turbulence is described. Some other interesting experimental results 
have been described by Lee and Ribner (ref. 9) which totally differ in principle from those 
of references 5, 6, and 7, in that the noise sources in a jet are described in terms of point 
measurements obtained by direct correlations between turbulence and sound; these corre- 
lations are not produced by the large-scale orderly structure but by a multi-eddy struc- 
ture of the flow. 


SYMBOLS 


A,B 

a 

b 

C 


c 

D 

T 

f 

f(t) 

G(r,t|r 0 ,t 0 ) 

h 

I(t) 

K 

k 

k l> k 2’^3>^4 ,k 5 


constants (eqs. (Bl)) 

longitudinal distance from jet axis to plane 

radial distance from jet axis to plane 

constant (eq. (8)) 

speed of sound 

nozzle diameter 

normal stress on surface 

frequency 

source strength 

Green's function (eq. (14)) 

source distance from the plane • 

acoustic energy flux 

pressure gradient transducer sensitivity factor 

wave number 

constants 
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L 

£ 

P 

P t 

Pm 

Po 

p(r,t) 


P^rms 

R 

R',R" 

R ff 

R V,Vp 



R Vp,Vp 


r 



s v,v P 


longitudinal distance from jet exit 
longitudinal distance from virtual origin (eq. (8)) 
acoustic power 
total power 

mean pressure 
static pressure 

fluctuating pressure 
root-mean-square pressure 
radial coordinate 

far-field radial coordinate (eq. (30)) 

correlation coefficient for point source strength (eqs. (A5) and (A7)) 

correlation coefficient of normal component of velocity and pressure 
gradient 

correlation coefficient for far-field sound 
correlation coefficient of pressure gradient 
position vector for point in space 
position vector for point on plane 
plane areas 

Fourier transform of I (eq. (20)) 

Fourier transform of Ry yp (eq. (22)) 
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S Vp,Vp 

S(x,t) 

Sff 

s 

rms 

^p 

d 

d' 


T t t 
U 3 

U i’ U j 

7 

V' 

i; 

V £ 

v n 

V r 

X, Y,Z 

X I , Y 1 ,Z 1 

X i’ X j 

x o 


Fourier transform of R Vp Vp 

transducer sensitivity 

Fourier transform of (eq. (A7)) 

root-mean-square pressure of pressure gradient transducer 
correlation function for far -field intensity 
far -field intensity 

far-field power per unit length of jet 
time 

jet velocity 

.instantaneous velocities 

1 i 

I 

mean square velocity 

1 

j velocity vector j 

longitudinal mean velocity (eq. (10)) 
velocity component normal to the plane 
radial mean velocity (eq. (9)) 

Cartesian coordinate system with respect to plane 
Cartesian coordinate system with respect to nozzle 
coordinate tensor system 

distance between virtual origin and nozzle exit 
angle between plane and jet axis (eqs. (37)) 



a 2 

0 

y 

6 

V 

e 

d 

H(a>) 

i 

Pm 
Po 
P( t) 

T 

$ 

’<b 

* 

U) 


w 


and 


angular relation of plane (eq. (36)) 
far -field angle with respect to nozzle 
interniittency factor 
delta function (eq. (14)) 
ratio, r/C 

far -field angle with respect to plane 
far -field angle with respect to plane 
phase angle 

distance between two points in the plane 
mean density 
static density 

fluctuating density ' 

time delay 

velocity potential (eqs. (5) and (8)) 
angle of plane 
mean flow angle (eq. (11)) 
angular frequency, 2rf 
dummy variable 

Primes denote derivatives with respect to variables. Asterisk (*) denotes conjugate 
< y denotes ensemble average. 
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THE FLUID FLOW PROBLEM 


Analysis 

Analysis of fluctuating flow outside the jet to determine the position of geometrical 
measurement plane .- Consider a fluid contained within a half space R bounded by a 
plane S. When there are no sources of energy within R and no transfer of heat by way 
of conduction or radiation, the rate of change of energy of the fluid within R is given by 
the total instantaneous energy flux across the boundary S (that is, the work done on the 
fluid inside R by the fluid outside R through the action of the surface stresses on S): 

I(t) = f (f • v) dS 
J S 


where V and f are the fluid velocity and stress at the surface. If a perfect inviscid 
fluid is assumed, (f • v) dS = pV n dS, -where p is the pressure and V n is the compo- 
nent of the velocity normal to the plane. The total energy flux then reduces to 



If this expression is to represent only the acoustic energy flux normal to the plane S 
coming from the jet, the plane must be positioned in a region where the acoustic pres- 
sure and fluctuating velocities predominate the mean and turbulent pressures and 
velocities. To answer the question of how closely one can approach the jet and still 

! i / 


remain in the acoustic region, the limit of the region where the wave equation' (that is, 

2 ^ \ ^ 

V 2 p = -4- + o(l 0 - ^)j ) is reasonably well satisfied must be determined. 
c 2 at 2 1 


This problem can be formulated by considering the basic flow-field equations and 
observing the behavior of the field in regions near and far from the jet. The continuity 
and the momentum equations for an inviscid fluid are 





( 1 ) 


The velocity, pressure, and density are decomposed into static, mean, and fluctuating 
components: 
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(2) 


u i = U i (r)+ u^r.t) 

p = P 0 + + / 

P = P 0 + P m (r)+ p(r,t) 


By neglecting quantities of higher order 1 and averaging over time, the mean components 
can be eliminated to yield 


2 fi ( M ) + p !^ = 0 
at ax i 

Pu aui(7,t) , apfa) - 


at 


ax; 


(3) 


for the fluctuating components. It is clear from the second of these two equations that 
the fluctuating velocity u^(r,t) is ir rotational. Thus, the fluctuating field described by 
equations (3) can be called "the irrotational field.” 


This irrotational field is now divided into two regions: region (1) close to the jet 
(in terms of acoustic wavelengths) where the velocity fluctuations appropriately normal- 
ized by the speed of sound as uj(r,t)/a Q are much larger than the relative density 
fluctuations p(r,t)/p Q and thus the fluctuating flow field can be considered incompress-j 
ible, and region (2) which is further away from the jet where the similarly normalized 
velocity fluctuations and the density fluctuations 'are of the same order, that is, 
p(r t) u i(r,t) 

— . In region (1) the continuity equation for the fluctuating flow reduces to 


9u,(r,t) 

Pq ^ =0 


8Xj 


(4) 


By introducing the velocity potential u(r,t) = V$, Laplace's equation 
V 2 $ = 0 


(5) 


* Here the perturbations of the mean quantities Uj(r), p(r), and p m ( r) have 
been assumed to be small, as well as the fluctuating perturbations u^(F,t), p(F,t), and 
p(r,t). 
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is obtained. In region (2) the density and the pressure fluctuations are related through 
the isentropic relationship 


9p(r ,t) _ 2 8p(r ,t) 
at at 


( 6 ) 


Again by introducing the velocity potential, the continuity and momentum equations 
reduce to the wave equation 


a 2 $ 

at 2 


- c 2 v 2 $ = 0 


(7) 


The region in which the Laplace equation is satisfied has been investigated by 
O. M. Phillips (ref. 10) and discussed by R. W. Stewart (ref. 11). These references pre- 
dict a decay of the mean-square fluctuating velocity with the fourth power of the distance 
from the jet. Measurements by Kibens (ref. 12), Bradbury (ref. 13), and Bradshaw 
(ref. 14) confirm this prediction. In the region farther away from the jet where the wave 
equation is applicable, it is well-known that the mean-square velocity decays with the 
second power of the distance. Hence, by measuring the decay of the mean-square fluc- 
tuating velocity, one should be able to determine an optimum position for the plane S at 
as small as possible distance from the jet boundary. Outside S, of course, the wave 
equation is then reasonably well satisfied. 


Mean flow quantities outside the jet .- By using the same arguments as in the pre- 
ceding section, it is readily seen that the mean flow in the irrotational region surrounding 
the jet satisfies Laplace's equation. Stewart (ref. 11) obtained a solution for the velocity 
potential of this axisymmetric conical mean flow which can be represented as 


$ = -C log e 


(( 2 + r 2 ) 1/2 - t] 


(8) 


where C is a positive constant and £ and r are the longitudinal and radial distances 
as measured from the virtual origin of the jet. The radial and longitudinal mean velocity 
components, V r and Vg, and the mean flow angle i// are then 


_ 3 $ _ 

.Cl 

1 + 1 

r ar 

I r\ 




(l + t). 

_ 94 > 

c 

1 

“ dSL ~ 

f (l ♦ , 2 ) 172 



(9) 

( 10 ) 
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where 77 is the ratio of the radial to the longitudinal distance r/£. 

Experimental Setup and Calibration 

The position of the plane used to evaluate the acoustic energy density flux along. the 
jet axis was established from a survey of the flow field about the jet. Measurements 
were made to determine the variation of the mean flow angle and of the mean and fluc- 
tuating velocities with distance from the jet boundary. 

In order to approximate the ideal situation of a jet expanding into an unbounded 
and otherwise quiescent region with no mean gradients, the test was set up in a room 
10 meters by 10 meters by 30 meters (fig. 1), that is, large enough to insure that the 
induced flow streamlines were not affected by the walls. The room acted as a huge set- 
tling chamber for the induced flow, which passed through a pair of wire damping screens 
in the upstream section. of the room. The nozzle consisted of a long pipe with a diam- . 
eter of 0.0625 meter connected upstream with an acoustic muffler that acted as a settling 
chamber for the flow. The flow from the jet was cold and had a mean exit velocity of 
Uj = 213 m/sec. ... 

All the flow measurements were made by use of hot-wire anemometry. Since many 
of the velocity measurements were extremely low (of the order of 1 m/sec), it was nec- 
essary to develop a new procedure for calibrating the hot wire at these low velocities. 

The rig used, shown in figure 2, consisted of a tube containing a propeller mounted on an 
essentially frictionless jewel bearing. 

By passing known quantities of air through the tube, the frequency of rotation of 
the propeller was related to the velocity of the air. Then the hot wire was placed in the 
tube and its output was calibrated against the output of the propeller. To insure relia- 
bility of the calibration, the upper velocity range was checked with standard calibration 
devices at the NASA Langley Research Center and at the National Bureau of Standards. 

Flow Measurements 

Inter mittency .- In order to determine the spreading angle of the flow, the location 
of the boundary of the jet was established from measurements of the intermittency y 
across the jet. The output of the conventional hot-wire anemometers placed in the flow 
shows irregular bursting of turbulence on a time scale as indicated by figure 3(a). A 
turbulence detector was used to discriminate between turbulent and nonturbulent flow. 

The signal was used to give the second time derivative of the velocity 1 3^v(t)/ 8 t^| ; this 
value converts into a square wave (fig. 3(b)) and then yields the value of y, that is, the 
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Figure 1.- Jet noise flow facility. 
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Figure 2.- Low- velocity hot-wire anemometer calibrator. 



Figure 3*- Signal conditioning with hot-wire anemometer. 







percentage of time that the signal is turbulent. As a result of the large-scale random 
convolutions of the jet boundary, the value of y varies smoothly from unity well inside 
the jet to zero outside the jet. 

With the boundary of the jet defined as the position where a particular value of 
intermittency occurs, the 1 -percent and 10-percent boundaries are shown in figure 4. 



Figure 4.- Jet spreading angle in terms of intermittency distribution. 


The figure indicates that the jet spreads at an approximately constant angle, namely, 11° 
for 1 -percent turbulence and 9° for 10-percent turbulence. A virtual origin of the jet, 
defined as the intersection of the boundary line with the jet center line, is seen to be about 
2.5 nozzle diameters upstream of the nozzle exit. 

Flow angle .- Measurement of the induced flow angle i// was made with a pair of 
cross-wire anemometer probes with a length-diameter ratio of 400:1. Figure 5 shows 
the results of these measurements at three stations along with intermittency measure- 
ments at the furthermost downstream station. It is seen that the flow angle gradually 
increases as the boundary of the jet is approached from the outside, reaches a maximum 
in the neighborhood of the boundary, and then quickly decreases to adjust to the direction 
of the jet. 
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0 . 1 1 10 

V 


Figure 5.- Flow angle and intermittency. t) = r/t. 

At large distances from the nozzle exit, the maximum flow angle approaches 90°. 
For this condition the surrounding fluid is entrained into the jet perpendicularly to the jet 
axis, which implies (ref. 2) that the rate of entrainment remains essentially constant as 
the jet progresses further downstream. This result conforms with previous measure- 
ments (Maestrello and McDaid, ref. 3) and indicates that the entrainment rate increases 
as the jet moves downstream until it reaches a constant value. 

Measurements of the mean flow angle at the furthermost downstream station agree 
well with those predicted by equation (11), whereas the upstream measurements indicate 
a smaller angle. (See fig. 5.) This condition is most likely to be caused by the physical 
boundaries of the flow. The analytical solution is based on flow issuing from a point in 
unbounded space whereas the experiment has the nozzle and upstream settling chamber 
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interfering with the induced flow. These physical boundaries have less effect further 
downstream where the analysis and measurements agree. 

Mean velocities .- Measurements of the mean radial V r and mean longitudinal Vf 
velocity components are shown in figure 6. These measurements are seen to follow 
equations (9) and (10) at the furthermost downstream station where the constant C is 
chosen to be 0.43 m2/sec. Measurements at the other two stations differ from the anal- 
ysis for the same reasons as the mean flow angle. The variations in the mean velocities 
at the three stations are similar; most noticeable is the decay of the radial velocity 1/rj 
with the inverse of the radial distance. The large increase in Vj at small tj indi- 
cates that these measurements were taken inside the jet boundary. This component will 
reach a maximum at the jet center line, whereas the radial component V r will fall off 
to zero. 

Fluctuating velocities .- The measurements of the mean-square fluctuating velocity 

u l^ + u 2^ + u 3^ are S iven * n fig ur e 7. As expected from the analysis, two distinct 
regions are noticed. The first is the incompressible region where the velocity fluctua- 
tions obey the Laplace equation and the mean -square velocity decays inversely with the 
fourth power of the distance. Farther from the jet center line, the mean-square fluctu- 
ating velocity decays with the second power of the distance, and thus indicates where the 
wave equation is applicable. The intersection of these two regions forms a boundary, 
which can be computed from the figure, and makes an angle of approximately 15° with the 
jet axis. These results provide information on how close the jet can be approached 
while still remaining in the acoustic field of the jet. In addition, the influence of the 
entrainment velocity on the acoustic energy flux will not pose a problem because the 
mean flow velocity is small and nearly uniform along the plane where the pressure is 
measured. 


THE ACOUSTIC PROBLEM 
Analysis 

The problem at hand is to express the acoustic energy density flux and the far- 
field acoustic intensity in terms of the near-field pressure gradient evaluated on an infi- 
nite plane at a distance as small as possible from the jet boundary, where the wave equa- 
tion is reasonably well satisfied. The minimum distance where the acoustic measure- 
ments should be made has been established in the preceding sections. The physical 
quantity to be used in the calculation is the gradient of the pressure, in particular, the 
measured correlation of this pressure gradient. 
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Fluctuating velocities in irrotational field 







on the plane.- The aim is to determine the acoustic energy flux normal to the element of 
area dS of a plane S inclined at an angle a to the jet axis (fig. 8) in terms of the 
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correlation of the pressure gradient on the plane. The acoustic pressure at any point on 
or above the plane (z ^ 0) (Morse and Feshbach, ref. 15) satisfies the wave equation and 
can be expressed as 




/- . i— 


dS o G ( r > t l r o’ t o) 


3p(r 0 ,t) 


3z 


(12) 


where the Green's function G satisfies the conditions 




The solution for the Green's function is known to be 


G ( r . t Fo» t o) 



1 

7tR 



du>’e 


• i R 



(14) 


where 6 is the Dirac delta function and R = r - r 0 |. Hence the acoustic pressure on 
and above the plane is given by 


^-oo R «J_oo 


-icu' ^ -iw'tj 
'e e 


9 P( r o> t+t l) 


9z 


(15) 


where 


t 


1 



- t 


2 

A generalized complex acoustic intensity in the n direction whose real part repre- 
sents, when r = 0, the normal component of the acoustic energy flux can be defined as 


l(7,r) = <p(r,t) V n *(r, t+r)> 


(16) 


By substituting the pressure from equation (15) into equation (16), the intensity becomes 


since 




dcu’e 


-iu>' — 
c 


-icu'tj 

e 




(17) 



9 P (Fq, t+4) \ 

9z / 



(18) 


By introducing the correlation coefficient between velocity and the pressure gradient, and 
on the basis of stationarity, one can show that 


2 

This generalized intensity which evidently is a correlation function of pressure 
with velocity is convenient when the subsequent equations are to be expressed in terms 
of correlations. 
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R v,v P (^» ‘<( v z‘(r’ t+T )|ffo' l )> * Tf ll VP’ t+r )|f( ;r 0’ t ) dt 


Therefore 


^ V z + t l) + ( T " t l)] ^ ° z _ R V,Vp( r ’ r o ; T_t l) (19) 


Taking the Fourier transform of equation (17) yields 


S T (r;w) = f l(r,r)e' iu;T dT 
J -oo 


i ^0 + n dS r>°° n°° • -iu>' — -iw't, 

= i? L dt i I s t?L d “' L dTe ' 1 “ Te Ce t -h) 


(20) 


Upon changing the variables of integration with the transformation 


t' = U 


t" = r - t, 


equation (19) becomes 


S .( 7 '“) ■ & I s TtII dt ' e '‘“ t ' II d “' e ’ 1 “ (c + ‘ ) II dt ” d ' 1 “‘" R V,Vp(?.?o;‘") 


where 


dS„ n0 + 


- I s -r L d ‘'e' l “ t ’M!+ ‘KvpP-V”) 

^I d ^"%vp(-» 


s v,v P ( 7 ^o'") * I__ dt e ' 1 “ tR v,vp( ? ’ r 0 ; t ) 


( 21 ) 


( 22 ) 


The experiments as discussed earlier yield the correlation of pressure gradient 
with a pressure gradient rather than its correlation with velocity. It can be shown with 
the use of the equation of motion that these two correlations are related by 
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r v,v p ( 7 ’V t ) = ^vp.vpfr^o; 7 ) 


8V z *(r,t) 3p(r o> t+r 
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az 


= lim -L f 

t^.oo 2T J. 1 


i r T 8V z ( F>t ) ap ( r o’ t+T ) 


at 


az 


dt 




= ifi 

x-»<x> 2T J_pp at 

1 pT *^v8p(r , t+T) 

fftnL v « 


V z *(7,t) 9 2 p(r 0 , t+r) 


at az 


dt 


_a_ 

ar 


dt 


" ~dT R V,Vp( r ’ r 0 ;T ) 


(23a) 


since 


where 


8V*(r,t) 3 p(7q. *t) \ t / ap »( 7 t ) 9p(f 0 , hr) 
at az / p V az az 


p R v P ,Vp( r ’ r o ;T ) “ “ dr R V, Vp( r,r o’ T ) 


K vp,v P V r » r o’ T -' \ az az / 




(23b) 




(24) 


By taking the Fourier transform of equations (23a) and integrating the left-hand side, one 
obtains 


S V,Vp( r ’ r o ;a; ) = ii S Vp,Vp( r ’ r o ;W ) 


(25) 


where 


S Vp,Vp( r ’ r o ;w ) " J ^ dre 1WTR vp,vp( r ’ r o ;T ) 
Hence equation (21) becomes 

‘“I 

s i (F;i,) * j s ds o *— s v P , vp( F - r oi") 


(26) 


(27) 
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Since the acoustic energy density flux through the plane is the main concern, r 
will be a point on the plane which, for convenience, may be chosen to be the origin of the 
coordinate system. In terms of polar coordinates £,$ in the plane, equation (27) 
becomes 


SjfOjw) 


1 

27 ?iu>p 



10J 


d£ e S Vp>Vp «U,<fco>) 


(28) 


where 


S Vp,Vp( 0> ^^’^ - S Vp,Vp(®’ r o ;a; ) 


The total power per unit length of the jet crossing the plane S is obtained by an inte- 
gration of the area elements dS around the jet. Since the jet is axially symmetric, 



where b is the radial distance from the jet axis to the plane. To show that equation (28) 
yields a correct analytical result, consider for the case of a monopole of strength f(t) 
located directly below the origin at z = -h. This case provides a numerical calibration 
of the method. The solution for the Fourier transform of the acoustic intensity from 
equation (28) is given in appendix A as equation (A8). 


S^Ojco) 



-i2( { 2 + h 2) 1/2 

e i|[(|2 + h2) 1/2 t {J 

(S 2 + l 
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This solution should equal the answer obtained by direct means (eq. (A9)) 


s l(0;u>) = (l + i 


C \ s ff M 
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Far-field intensity from correlation on the plane . - The far-field intensity will now 
be determined from measurements of the cross spectral density S Vp vpl^’^o’^) on ttie 
plane S. Since there are no sources in the volume above the plane, the far-field pres- 
sure is given by equation (12) which can be written as 


23 



where R = | r - r Q |, r Q is on the plane, and r is the far-field point as shown in 
figure 9. The cross correlation of the pressure at r" is 


Rp p(? ;T ) = <( p(r.t) p'(r, t+x)^ 


-i ttR' i^R" 


1 ri + +T /> n p°° p°° icWt-t,/) 1 1 / \ 

■nzL“° , L a vU«'£ a vL'‘"'L a "''* * 


R fi-T- t-M 8 p >~o"V) 

n Vp,Vpl r o ,r o » L o ~ l o) ■ dz dz 



Figure 9*- Far-field geometry for calculation of the intensity. 
dA = (r sin d d0) (r dd) . 
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Without any further loss of generality, one might let t = 0 and introduce t Q = t " - r; 
then the spectral density of the mean-square pressure at r is 


dTe' 1WT R pp (?;T) 


-L 

D - v is '*■£>£> 




■ R vp,vp( ? o'- r o"; T+t o-‘o') 


Changing the variables of integration by transforming 


(31) 
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yields 


s n „(?;«) - f dt, f dt, f dt. f dS ' R_ f dw'e ' 2 c ' f 

PP 16.4 3J_„ 2J_„ ljo o R'R" Vp.VpV o > o - 3; 


iu;' 


du>" e 




, r°° rO* r0 + r n " iwt 3 iwt l ~ iut 2 , 

= ^2 L dt 3 L dt 2 L dt l Is “o’ Is dS o" 5 Wr^ R Vp,Vp( F o’> r o"- t 3) 9 ( t 2 + “Hh + T~) 


which reduces to 
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(32) 


where 


s vp,v P ( ? o’’ r o , '; u ’) = dTe " 1 " TR vp,vp( i: 'o'- r 0 "; T ) 


Now consider r to be the far-field point shown in figure 9 and express it in terms of 
the spherical coordinates r , 3 , 9 . Thus, 


R' = r 



sin 3 ( x , CQS q + y, sin + o 
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and 


lim 

r-*oo 


i|R' -ifR" 

.2 e c e c 


R' 


R” 


i^sin aJ(x ,, -x , )cos 0+(y"-y')sin ej 


As a result, the spectral density of the mean-square pressure at the far-field point r 
becomes 


oo oo 

s P,P (r - e - 3 :“) = 13 ff dx ’ dy ’ If dx " dy " W^") 

_oo _ oo 


i^sin a[(x’'-x')cos 0+(y”-y')sin 9 


(33) 


One can calculate the contribution to the far-field spectral density of the pressure 
coming from the point r Q which, for convenience, is chosen to be the origin. Define 

oo 

S p,p( r,0 ’^ ;w ) = JT dx ' 4y'^(r,M.r 0 , ;w) (34) 

- oo 

By using the notation of figure 8, that is, 

x" - x' = 4 COS <p 
y" - y' = I Sin 0 


one obtains from equation (33) 

(2 ttt) 2 dir, 6,3,0 ;w) = f d<f> f £ d£ e 


i -±- 1 sin 3 cos(cp-d) 


s vp,vp(°’Z’to u ) 


(35) 


In effect, the far-field spectral density has been replaced by a sum of contributions from 
apparent sources at each point on the plane S. It is of more interest to calculate the far 
field in term of contributions per unit length of the jet. Since the correlation distribution 
is the same at each circumferential point of an axisymmetric jet, the plane appears to be 
rotated about the jet to form a conical surface wrapped around the jet. By referring to 
the notation of figure 10, the contribution to the spectral density of the mean-square 
pressure at the far-field point R,/3 coming from the conical strip at x^ = a can be 
formulated as 


^'(R, ft to, Ax) 


I A r > 

Strip 


0,3, 0;co) dag 


(36) 
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Figure 10.- Far-field geometry used for determining the power per unit length. 

It now remains to express the far-field coordinates with respect to the plane r,0,3 in 
terms of the far-field coordinates with respect to the jet R,/3 and the angle a 2 at 
which the plane is rotated. The Cartesian transformation between the two far-field coor- 
dinate systems is obtained by simple translation and rotations as 

x = (xj - a)cos + (z^ cos a 2 - y^ sin a 2 - b)sin a 

y = y^ cos q »2 + b sin ce 2 cos a 2 + z i s ^ n a 2 cos a 2 


z = -^Xj - ajsin a + cos <x 2 " sin a 2 ~ b)cos a 


By noting that for the far-field point of interest, 


Xj = R cos f3 
y 2 = -R sin /3 
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one obtains 


x = R(cos p cos a + sin p sin a sin a^) “ a cos a ~ b sin a 
y = -R(sin p cos a^j 

z = R^sin p cos a sin - cos p sin a) + a sin a - b cos a 


Letting R become large results in 

d/2 


r = (x 2 + y 2 + z 2 ) 1 ^ = R 

0 = cos“*(z/r) = cos -1 (sin p cos a sin ffg - cos p sin a) 




(37) 


6 = tan' 


'‘( I )-*”' 1 


-sin p cos d 2 


cos p cos a? + sin p sin a sin a 2 


j 


The limit of integration in equation (36) is obtained by noting that the plane S will be 
below the far -field point only when z > 0. This condition results in 


sin dig 


> tan a 
tan p 


Hence, the contribution to the spectral density of the pressure at the far-field angle p 
with respect to the jet axis coming from the strip at x = a is given by 

* 

p7T-Q?2 / — v 

^’(R.frw.AXj) = \ + ^(r,0,a,O;wj da 2 

a 2 

where r, 9, and 3 are given by equations (37) and 



System Calibration 

In order to validate the experimental procedures in the measurement of the acoustic 
intensity, a test was conducted with an experimental point source. The design of the 
point source was similar to that of the one developed by Atvars, Schubert, Grande, and 
Ribner (ref. 16). The experiment consisted of evaluating Sj(0;u>) (eq. (28)) from 
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measurement of the correlation function S Vp vp (o,£,0;u>) from the experimental point 
source. This correlation function was measured by using two pressure gradient micro- 
phones. The development and calibration of the microphone is presented in appendix B. 
The result was then compared with the measurement made with a single condenser 
microphone system placed in the vicinity of the source which determines the intensity 
directly. The cross spectral density from a point source along (p = 45° at a frequency 
f = 1809 Hz is shown in figure 11, and the extreme end of the spatial correlation extends 
up to R e = 50 cm. 

The intensity becomes 


s i (0: "> * zirap J" e dR i £ /2 d,p s vp,vp(°> R i'*;") cos (^) iS v P ,v P (°' R i'*;") 


The correlation function of Sy p y p as obtained from the curve fitting of figure 12 is 
represented by 


, “ k l^l 

S Vp,Vp = k 5 e cos 


(ko^ sin <p + ko^ cos ^ / o\ 

“ 


and 


kj = 0.051 
k 2 = 17.0 
k 3 = 13.20 
k 4 = 0.053 
k 5 = 1.890 

The comparison between the results obtained from the integration of the measured corre- 
lation with the direct reading of the intensity was within 22 percent, or about 1 dB. This 
difference is the order of magnitude of the cumulative error to be expected between the 
two methods used to evaluate the intensity when all possible errors of the two systems 
of instrumentation which include transducers, auxiliary equipment, correlator, and read- 
out equipment are estimated. 
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STATION 5 
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Figure 12.- Broad-band correlation of pressure gradient. 


Acoustic Measurements 


Correlation measurements were made on a geometrical plane placed at a distance 
31.5 cm from the center of the nozzle, at an angle a = 17 o 20' relative to the center 
line of the jet (fig. 13). The distance from the plane to the jet was chosen from the mea- 
surement made of the velocity field so that the pressure on the plane is an acoustic 
pressure. Measurements of the space-time correlation of the pressure gradient 
R Vp V p(0, were made within the band between 250 Hz and 3000 Hz at 10 fixed 

points along the plane at L/D ratios between 0 and 20.56 lying on the projection of the 
center line of the jet. The correlations were made along three axes lying on the plane 
and passing through r . The axis corresponds to the rays of 0=0 and ir; tt/A and 
577/4; and 0 = 7 t/ 2 and 37 t/ 2 in figure 8. Correlation measurements were made for 
10 stations along the plane. Typical correlation measurements along these three axes 
are shown in figures 12, 14, and 15 for station 5 and correspond to L/D = 9.26. For 
higher frequencies, these measurements allowed R Vp Vp (0, 4,0 ;t) to be defined with 
sufficient accuracy whereas for lower frequencies, additional test points on a greater 
number of axes would be desirable to improve resolution. It is clearly evident that the 
sound field is nonhomogeneous; this statement implies that as the sound field is radiated 
away from the jet, it has a preferred downstream direction. The cross correlation is 
the same along the rays 0 = 0,77 and 0 = 77,0 because the flow from the jet is also 
axisymmetric; therefore only the correlation along 0=0 is shown. The spatial extent 
of the cross correlation indicates that the unit area of the plane is sensitive to fairly 
extended regions of the jet especially at large L/D values where the lower frequency 
predominates the higher. For the stations downstream of station 5, the spatial extent of 
the correlation R Vp Vp along the rays 9 = 77/4 and 77/2 decreases less rapidly with 
distance than along the rays 9 = 577/4 and 377/2 whereas for stations upstream of sta- 
tion 5, the correlation R Vp Vp becomes more homogeneous in space. These results 
imply that for the stations close to the jet exit, the sound field radiates away from the jet 
nearly omnidirectionally; whereas for stations further downstream, the sound field is 
directional. This picture of the correlation field describing the acoustic wave pattern as 
it emanates from the jet is consistent with the far -field intensity measurements in this 
experiment and with results of other experiments as described in the literature. 

The cross spectral density S Vp y p (0, £, 0;o>) was obtained by taking the Fourier 
transform of the broad-band cross correlation for 12 center frequencies between 250 Hz 
and 3000 Hz at 250-Hz intervals. From the narrow-band correlations, one can construct 
the wave vector spectrum as was done previously by Maestrello and McDaid (ref. 2). 

This spectrum displays a pronounced maximum on wave number space which corresponds 
to the direction of the sound reaching the plane of measurements. Although this spectrum 
will not be reported, it can be obtained from the available information. 
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STATION 5 



Broad-band correlation of pressure gradient. 




STATION 5 
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Figure 15.- Broad-band correlation of pressure gradient. 





An additional measurement of the pressure field was made at one fixed location 
along the jet by using two condenser pressure microphones to insure that the choice made 
on the position of the plane along the jet is consistent with the measurement made of the 
velocity field. This measurement shows that the velocity along the plane satisfies the 
acoustic wave equation. The result shows that all the convected components except that 
for the frequency 250 Hz of the cross correlation about the geometrical plane are acoustic 
components . 


PHOTOGRAPHIC REPRESENTATION OF SOUND FIELD 
PATTERN EMANATING FROM A JET 

Noise radiation from jet exhausts can be represented by a pattern of colors, 
obtained by means of a photographic plate, each color corresponding to a different inten- 
sity level of noise. This effect is seen in figure 16 which represents the noise field from 
a jet in that region where the wave equation is reasonably well satisfied. Although this 
is a new application in aeroacoustics, the method is analogous to a color image formed 
in a television picture tube. This photographic technique was first developed by Kock 
and Harvey (ref. 17) and by J. T. Rainey and D. G. Neville (ref. 18). The colors in the 
picture are the result of time exposures of five different color light bulbs, each bulb being 
triggered by the output voltage of a nearby condenser microphone mounted on a swinging 
arm. The microphone traverses successive arcs arranged in such a manner that their 
loci are a vertical plane at an angle 17°20' to the thrust axis and at a distance of 31.5 cm 
from a 6.25-cm-diameter nozzle. The resulting figure 16 is made up of superimposed 
elliptically shaded color patterns stretched along the flow axis. The maximum intensity 
region is represented by red; and orange, green, and yellow represent successive 2-dB 
decrements in noise level. 

This technique should be developed further because the display of the sound field 
patterns is a useful representation of the distribution of the sound pressure level. This 
technique can easily be extended to represent sound energy, energy density flux as well as 
propagation path, and all quantities required to evaluate the mechanism of sound genera- 
tion and propagation. 


RESULTS AND DISCUSSION 

Evaluation of the Acoustic Energy Density Flux 

In order to insure the accuracy of the techniques employed in this analysis, com- 
parison is made in this discussion between the total acoustic power and the acoustic 
intensity, evaluated by equations (28) and (36), respectively, with the corresponding quan- 
tities directly measured by a simple microphone in the far field. New information of 
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L/D = 54.4 



Figure 16. - Broad-band jet noise field visualization. U. = 213 m/sec; lens width, 166 cm. 
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acceptable accuracy about noise emanating from a jet is presented, namely, the distribu- 
tion of acoustic power along the plane of measurements and the far -field intensity result- 
ing from the contribution of sound emanated from the various stations along the plane. 

Evaluation of the acoustic energy density flux Sj(0,f) (eq. (28)) is shown in fig- 
ure 17, for 10 stations along the plane. The flux is normalized with the mean value 

/-»QO 

( Sj(0,f) df at each station and is plotted in terms of the Strouhal number fD/Uj and 

in terms of L/D, where D is the jet diameter, Uj is the jet exit velocity, and L is 
the distance from nozzle exit to the projection of the plane into the jet axis. (See fig. 13.) 

It is evident that the higher frequency components are contributed mostly from the region 
near the jet exit and this contribution diminishes rapidly with increasing distance down- 
stream, whereas lower frequency components, mainly 750 Hz and 1000 Hz (corresponding 
to fD/Uj = 0.225 and 0.3), are contributed more evenly from all stations along the plane. 
These lower frequencies also correspond to the peak amplitude power spectral density in 
the far field. 

Some caution must be exercised in interpreting the source distribution from the 
evaluation of the acoustic energy flux. The distribution of the acoustic energy density flux 
on the plane near the jet cannot be interpreted as an equivalent source distribution. What 
is evaluated on the plane is relevant only to the sound field as it emanates from the outer 
boundary of the jet and not in the manner in which the source originates sound since the 
source inside the jet has finite life time, and the distribution of the radiation depends on 
the resultant instantaneous strength by the superposition of number of multipole fields. 

As the pressure is released from the source, it interacts with the jet stream extracting 
energy therefrom as well as imparting energy thereto. The full convected wave equation 
provides the necessary clue as to how the pressure is transmitted through the flow. The 
sources in the jet appear and disappear as they move imbedded in the moving stream and 
indicate that in this mechanism, the steady-state condition does not exist . This is one of 
the main reasons why the sound field is evaluated outside the jet before it is evaluated from 
the inside. 

The broad-band emission of sound per unit length of the jet evaluated on the plane is 
shown in figure 18. It is obtained by rotating the plane around the jet and adding the con- 
tribution made per unit area of the plane. Since the jet is axisymmetric, it is assumed that 
the power per unit area of the plane is the same at the same radial position of the plane 
with respect to the jet. The sound emitted peaks at approximately 9 jet diameters down- 
stream of the jet and it varies in amplitude by a factor of better than two at the extreme 
stations. Intuitively, one would expect the distribution of sources inside the jet to have a 
distribution which is shifted more closely to the jet exit than the one measured on the plane 
due to convection and refraction as the sound crosses the flow on its way out of the jet. 
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Figure 17. - Normalized acoustic power per unit area of the plane 
between stations 0 and 9- 
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In order to insure that the area inside the curve of figure 18 corresponds to the 
correct value of the total acoustic power, comparison was made with far-field total power 
measurements between 250 Hz and 3000 Hz. The direct far-field measurement was 
obtained by integrating elemental areas around a 3.65 -meter-diameter sphere. The 
results show that the total acoustic power from far-field measurements is 0.82 watt and 
the total acoustic power from figure 18 is 0.74 watt. 

The comparison is remarkably close about 0.5 dB between the two approaches; it 
validates, at least in part, that this technique is of practical significance in obtaining the 
distribution of the acoustic power along the jet. Additional information on the distribution 
of the acoustic power along the jet can be obtained by plotting the power per unit length 
(normalized with the total jet power) against Strouhal number fD/Uj (fig. 19). Peak 
amplitudes occur for fD/Uj = 0.3 and 0.225 which also is the range of Strouhal numbers 
where the far-field power spectral density peaks. The variation of the power with dis- 
tance is very gradual except for the initial few diameters. For much higher Strouhal 
numbers the distribution shifts close to the jet exit and the field is not as broadly spaced 
as that for low Strouhal numbers. Therefore, the significant medium frequency far-field 
contribution comes from a region of the jet conical outer boundary where the emission of 
sound is stretched well beyond 10 nozzle diameters and is not as localized as those for 
the higher Strouhal number components. However, as none of the components show strong 
localized regions of sound emanating from a particular part of the jet, this result tends 
to suggest that there are no highly localized directional sources in the flow. These 
results relate the distribution of sound emission on the plane and ought not to be com- 
pared with the distribution of sources obtained by using similarity considerations by 
Ribner (ref. 19), Powell (ref. 20), Lilley (ref. 21), nor with recent measurements of direct 
correlation made by Lee and Ribner (ref. 9) since one refers to the source inside the jet, 
whereas this work refers to the distribution of sound as it exited from the jet. 

Evaluation of Far -Field Intensity 

Evaluation of far-field intensity A'( R,/3;o>,Ax) (eq. (36)) was checked out by using 
a known function for Syp yp(o ,£, 0;o>) with omnidirectional spatial distribution over the 
axis of the plane. The far-field intensity A' shows an equivalent omnidirectional inten- 
sity except in the region where the angle /3 approaches the angle of the plane a for 
reasons presented in the formulation of the problem. This comparison validates the 
numerical procedure used in evaluating A' ■ Computation of the function A' was made 
at two far-field angles, 90° and 45°, and the results are plotted in terms of the Strouhal 
number fD/Uj in figures 20 and 21. In addition, the far-field intensity obtained by add- 
ing the contributions to the intensity evaluated along each station up to L/D = 20.96 was 
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computed from equation (36) for different locations. 







Contribution to far-field intensity at 3 = 45 , R = 3-65 meters along the plane 
computed from equation (36) from different locations. 






compared with a single microphone measurement in the far field (fig. 22). The agree- 
ment for both 90° and 45° is satisfactory. 

The additional contribution by intensity in the region of the jet beyond L/D = 20.96 
is not significant, as deduced from figures 20 and 21. In this experiment the minimum 
acceptable far -field angle /3 is approximately 30° and is limited by the plane angle a. 
One can reduce the plane angle, however, by increasing the distance between the plane 
and the axis of the jet. (See fig. 13.) 

The distribution A' is the far-field intensity as made up of contributions from 
each station on the plane plotted in terms of the projection of the plane into the jet 
axis L/D. Conclusions can be drawn from figures 20 and 21 regarding what is the most 
significant region of the jet contributing to the far field for a given Strouhal number. It 
is significant from the figures that fD/Uj = 0.30 and 0.60 contributed mostly at p= 90°, 
whereas at /3 = 45° most contribution comes from fD/Uj = 0.30, 0.225, and 0.60, and 
indicates an apparent shift in predominant spectral distributions to the lower frequencies 
between 90° and 45°. 

The distribution d' for /3 = 90° has a more clearly defined maximum value as 
L/D varies than that for 0 = 45° which shows a more extended contribution along the 
various stations on the plane. The intensity amplitude between the two far-field angles 
varies in a systematic manner, decreasing between 90° and 45°. The directivity and the 
intensity at each station along the plane is dependent on the cross correlation of the pres- 
sure on the plane, and since the plane is relatively close to the jet boundary, these results 
provide an initial clue of the contribution of each station of the jet to some particular far- 
field angle. This result shows that the emissions of sound from various stations differ 
more with far-field angle than with Strouhal number, as far as the spatial extent of emis- 
sion at various stations is concerned; the Strouhal number only adjusts the relative 
amplitude of emission. This result may clarify some directivity anomalies reported by 
Lush (ref. 22). Another interesting observation is that there appears to be some signifi- 
cant shifting of the peak of the distribution with far-field angle and with Strouhal number; 
this trend, however, is not well characterized. It could be due to some inaccuracy of the 
data between adjacent stations; any inaccuracy, however, is not expected to alter dras- 
tically the observed behavior. This last observation may be more apparent since the 
change in far-field angle tends to compress or stretch the contribution made by each sta- 
tion to the far-field point. 


CONCLUDING REMARKS 

This paper establishes the relationship of sound emanating from various regions of 
a well-defined conical outer boundary of a subsonic jet with the far-field directivity and 
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intensity distribution. The technique used permits acoustic measurements to be taken 
close to the regions from which the sound emanates without getting involved in detailed 
measurements of the flow field inside the jet. 

To define the geometrical plane where acoustic measurements were carried out 
and which, in turn, define the conical outer boundary, velocities were obtained by mapping 
of the flow field about the jet outside the vortical flow. This mapping shows two distinct 
regions of interest: the region close to the interface, where Laplace's equation is valid, 
and the mean-square value of the velocity decreases inversely to the fourth power of dis- 
tance, as originally predicted by Phillips, and a region further from the jet, where the 
acoustic wave equation is satisfied with the mean-square velocity decreasing inversely to 
the second power of the distance. The induced flow approaches the jet at approximately 
90° from the flow axis and the mean radial velocity component decreases inversely with 
the first power of the distance. These results are consistent with the prediction made 
by Stewart. 

The results show that the broad-band acoustic power evaluated on the geometrical 
plane peaks at approximately 9 diameters downstream from the nozzle. Also, the 
narrow -band acoustic power at a Strouhal number of about 0.3 predominates over nearly 
the entire length of the jet. The contribution of the sound from the plane of measure- 
ments to the far field, however, depends less on the Strouhal number than on the field 
angle. This result may clarify some directivity anomalies reported by Lush. 

One can observe that the far-field intensity distributions along the jet differ signifi- 
cantly for the angles of 45° to 90°. At 45° the intensity distribution for Strouhal num- 
ber 0.9 is far below the mean level of intensity for Strouhal number of about 0.3, whereas 
at 90° for the Strouhal number at 0.9, the intensity distribution is much closer to the 
peak distribution. 

These measurements show that the far-field frequency spectrum peak is related to 
the spatial extent of the region where the sound emanates from the jet. The 45° com- 
ponents emanate from an extended region of the jet, over 17 diameters long, whereas the 
90° components of the sound field are concentrated in a narrower region. It is apparent 
that a large increase in jet velocity cannot significantly alter the spatial extent of the 
region from which the sound emanates at 45° unless the flow becomes supersonic. The 
-change will affect the amplitude distribution but will not shift the peak frequency since 
the amplitude of the higher frequency components cannot change enough relative to that 
at a Strouhal number of about 0.3 to alter the peak. However, for the 90° components, 
the higher velocity can easily alter the peak frequency since the contributions from the 
different Strouhal numbers are relatively close together. Thus, it may again be inferred 
from the present investigation that the far-field spectrum peak depends more on the far- 
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field angle than on the Strouhal number for far -field angles of the order of 45° or less. 
The Strouhal number dependence of the normalized spectrum distribution becomes sig- 
nificant only for large far -field angles. 

These results can also provide information relevant to sound suppression by means 
of surface shielding or shielding by refraction, obtained by tailoring the temperature gra- 
dient of a jet, in terms of the length of the surface or that of the refraction zone. This 
conclusion is drawn from the spatial extent of the region where the sound is emanated 
from the jet. At the 45° far-field angle, substantial suppression can only be achieved 
when the shielding surface or the refractive zone extends beyond 20 nozzle diameters 
downstream and beyond 10 nozzle diameters laterally as seen from the lateral cross 
correlation whereas for 90° far-field angle suppression can be achieved by using a 
smaller surface. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 23, 1973. 
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APPENDIX A - Continued 


By introducing the variables 
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By taking the Fourier transform, 
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where 


S ff M=f dTe" iwr R ff (T) 


Hence, from equation (28), 


^ C d4 d * ^ c ( s vp,vp(». 


r 

__ s _ r e ‘^ s tf M ft + ‘ I 1 -) 6 ' 1 '* 1 , [■ - + h2 ) 1/2 } 

iU} P •'0 ($ 2 + h 2 )_ h J|^ (^2 + h 2 ) 


1/2| i^(? 2 + h 2 ) 1/2 


. ft * *S») 


• U) L 

-l^-h 


e s ff(“> r” 

iwph J 0 * 


1 . i|({ 2 + h 2 ) 1/ je i ^^ 2+h2)1/2+ ^ 


(« 2 + h 2 )' 


The exact answer for the intensity is found by considering 


l(0;r) = <(p(0, t + r)v z + (o,t)^> 


Thus, differentiation of this expression gives 




=K p(5 ' 


9plO,U 

9z 


Substituting the pressure and gradient from equations (Al) and (A2) 


pch“ ' ' * 

1 11 r ( T ) 

pch2 9 t ph3 ff 


t, = t - * 
1 c 
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The Fourier transform of the intensity is 


Sj(0;<u)=f l(o ;t ) e" ia,T dr 

A v/ -OO 

= J_ f° 91(011) e -icur dT 
id* J.OO 9 t 



This equation should be equal to equation (A8). 


APPENDIX B 


DEVELOPMENT AND CALIBRATION OF THE PRESSURE 
GRADIENT TRANSDUCER 


The evaluation of the integrand S V p yp(o , £,$,u>) is made from measurements of 
the correlations of the pressure gradient in a plane near the boundary of the jet. The 
selection of a pressure gradient transducer over other types of pressure transducers 
was made because of its directional characteristic, its discrimination against distant 
sound sources, and its rejection of randomly incident background noise. The pressure 
gradient transducer, shown in figure 23, is of the capacitance type, an "electret" foil 
being used. It is similar to that developed by Sessler and West (ref. 23) with the excep- 
tion that it is smaller in size and has an adjustable diaphragm tension. This additional 
feature in the transducer design was introduced by Carl Rumble of NASA Langley 
Research Center. The diaphragm of the transducer, both sides of which are exposed to 
the pressure field, is a permanently polarized dielectret foil made of 0.0037- cm metalized 
Mylar. The case is made of brass and is about 1.25 cm in diameter and 0.96 cm thick. 

The calibration is made in an impedance tube for frequency response, sensitivity, direc- 
tivity, phase distortion, and harmonic distortion. The most reliable method of calibration 
is the one made in the impedance tube (fig. 24). Let the incident wave have an ampli- 
tude A and the reflected wave have an amplitude B. The pressure is given by 


p(x,t) = A cos(cut - kx) + B sinjurt + kx + p(co)] 

= |a + B cos[2kx + p(w)j| cos(wt - kx) - B sinjjlkx + p(w)]sin(u>t - kx) 


The root-mean-square pressure becomes 

[pW] - 1 


Irms 


(2)V2 


A 2 + B 2 + 


2AB cos(2kx + p(wj]^| 


(Bl) 


(B2) 


Then from equation (Bl) the pressure gradient becomes 


9 ^9x^ = ^ sin(cut - kx) - ^B sin[u>t + kx + /i(o>)J 

= ^ |a - B cos[2kx + p(u>)[|sin(ct)t - kx) 

- — B sin[2kx - p(cu)]cos(a>t - kx) 


(B3) 
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Figure 23.- Pressure gradient transducer. 
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Figure 2k. - Calibration of pressure gradient transducer. 



APPENDIX B - Continued 


and the root-mean-square pressure gradient can be written as 


9£M 


9X 


-Jrms (2) 
= 

- 


T/2 + 2AB cos|2kx + u (u))] - 4AB cos[2kx + u (w)J 


1/2 


rms 


rms 


'1 - 


4AB cos (2kx + /x(cu)] 


)l/2 


A 2 + B 2 + 2AB cos[2kx + p(a>j] 

(1 + e) 


(B4) 


where |e| g . 

A' 2 + B z 

From equation (B4) the error e made in the calibration of the pressure trans- 
ducer depends on the amplitude of the standing-wave ratio. To minimize this error, the 
downstream end of the impedance tube was treated with an anechoic termination. Calcu- 
lations show that for a standing-wave ratio of 0.25 dB, the amount of error introduced is 
of order e = 3 percent. The calibration was made for those frequencies in which the 
standing-wave ratio did not exceed this value. 

The sensitivity is defined as 

S(x,t) = K 
v ’ 9x 


and from equation (4) 


K = 


c 

rms c_ 
' 9p \ oo 
, 3x /rms 


(B5) 


The frequency-response calibration is made by placing the transducer in a spherical 
sound field generated by an artificial mouth. The response is shown in figure 25(a) where 
the curve is flat up to 3000 Hz. The transducer measured the component of the pressure 
gradient which is normal to the transducer face. This result is apparent in the figure 
where the sensitivity is seen to be down 3 dB when the gradient is incident at 45° to the 
face. The sensitivity factor K is plotted in figure 25(b); it remains constant up to a 
frequency of 3000 Hz. The phase-angle calibration was made by measuring the phase- 
angle difference between output voltage and incident pressure. The difference was 
within 3° from 200 to 3000 Hz. Furthermore, the harmonic distortion for an input pres- 
sure between 85 dB and 125 dB was linear within 5 percent. 
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Frequency, kHz 


(a) Transducer frequency response. 



Frequency, kHz 

(b) Transducer sensitivity factor, K. 

Figure 25.- Transducer calibration. 

Since this measurement was made, further developments have taken place in 
improving the frequency response and the sensitivity of the transducer. Additional 
studies pertinent to the design of the pressure transducer in terms of the mechanical, 
acoustical, and electrical characteristic have been recently completed by E. Gonzaga 
(ref. 24). 
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